Stripes and superconductivity in the two-dimensional self-consistent model 
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We found solutions of the Bogoliubov-de Gennes equations for the two-dimensional self-consistent 
model of superconductors with d x 2_ y i symmetry of the order parameter, taking into account spin 
and charge distributions. Analytical solutions for spin-charge density wave phases in the absence of 
the superconductivity ("stripe" and "checkerboard" structures) are presented. Analytical solutions 
for coexisting superconductivity and stripes are found. 

PACS numbers: PACS numbers: 71.10.Fd, 74.72.-h, 64.60.-i, 71.27. +a 



I. INTRODUCTION 



Stripe order, e.g. coupled spin- and charge- 
density periodic superstructure, is found in the under- 
doped superconducting cuprates and competes/coexists 
with superconducting order. Charge and spin 
stripe order have been observed experimentally in a 
some cuprate compounds, specifically La2- x Ba x CuO^ 
and Lax.&-xNdo.iSr x CuO^^. A recent study of 
Lai.s7bBao.\2bCuOi with angle-resolved photoemission 
and scanning tunneling spectroscopies^ has found ev- 
idence for a d-wave- ltype gap at low temperature, 
well within the stripe-ordered phase but above the bulk 
Tc. An earlier inelastic neutron scattering data 5 shows 
field-induced fluctuating magnetic order with space pe- 
riodicity 8ao and wave vector pointing along Cu-0 
bond direction in the ab-plane of the optimally doped 
Lai.84Sro.i6Cu04 in external magnetic field of 7.5 T be- 
low 10 K. The applied magnetic field (~ 2 — 7 T) imposes 
the vortex lattice and induces "checkerboard" local den- 
sity of electronic states (LDOS) seen in the STM exper- 
iments in high-T c superconductor E^S^Ca C^Os+s 6 - 
The pattern originating in the Abrikosov's vortex cores 
has 4ao periodicity, is oriented along Cu-0 bonds, and 
has decay length ~ 30 angstroms reaching well outside 
the vortex core. The existence of antifcrromagnetic spin 
fluctuations well outside the vortex cores is also discov- 
ered by NMR 7 in superconducting YBCO in a 13 T ex- 
ternal magnetic field. Theoretical predictions had also 
been made of the magnetic field induced coexistence of 
antiferromagnetic ordering phenomena and superconduc- 
tivity in high-T c cuprates^ - — due to assumed proxim- 
ity of pure superconducting state to a phase with co- 
existing superconductivity and spin density wave order. 
In these works effective Ginzburg-Landau theories of cou- 
pled superconducting-, spin- and charge-order fields were 
used. Alternatively, the fermionic quasi-particle weak- 
coupling approaches were focused on the theoretical pre- 



dictions arising from the model of BCS superconductor 
with d x 2_ y 2 symmetry^. An effect of the nodal fermions 
on the zero bias conductance peak in tunneling studies 
was predicted. However STM experiments of the vortices 
in high-T c compounds revealed a very different struc- 
ture of LDOS 6 . In this paper we make an effort to com- 
bine both theoretical approaches and present analytical 
mean- field solutions of coexisting spin-, charge- and su- 
perconducting orders derived form microscopic Hubbard 
model in the weak-coupling approximation. The previ- 
ous analytical results obtained in the quasi ID cases^— 
are now extended for two real space dimensions. Dif- 
ferent analytical solutions for collinear and checkerboard 
stripe-phases, as well as for spin-charge density modu- 
lation inside Abrikosov's vortex core are obtained. Si- 
multaneously, our theory provides wave-functions of the 
fermionic states in all considered cases. 



II. EFFECTIVE HAMITONIAN. 
BOGOLIUBOV-DE GENNES EQUATIONS 

Consider the Hamiltonian H = Hq + H sc consisting 
of two parts: the first part is the Hubbard Hamiltonian 
with on-site repulsion U > 



(i,j>,<7 i i-<? 

and the interaction part including superconducting cor- 
relations 



h sc = A ^j^) c lA-- + h -' 



(2) 



<i,j>,a 



where 2<j j> a ^ s a summation over nearest neighboring 
sites rj, rj of the square lattice, and spin components 
a = 2s, = ±1. 
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In the self-consistent approximation the Hamiltonian 
acquires the form 



H 



-E 

2 ^ 
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-UJ2(S,M)W4,^ + U{S z {r t )) 2 - ^E C L C - 



<i,]>,a 
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-,(3) 



where we introduce similar— slowly varying functions for 
spin order parameter m(rj) and the charge density p(?i) 
defined as 



p(r) = (n(r)), (-If* +«m(r i ) = U(S z (r z )), 

A(i,j;o-) = -g{c-j,-<jCi,a)- 



(4) 



We can diagonalize the total Hamiltonian H — Hq + H sc 
by performing a unitary Bogoliubov transformation 

C<r(r) = ^ 7n,aWn,a(r) - O^-o^-aO) ( 5 ) 



New operators 7, 7 + satisfy the fermionic commutative 
relations {7ri,er, 7„ o-/} = <5m,ri<5<T,er'- The transformations 
([5]) must diagonalize the Hamiltonian H : 



(6) 



when (— \~) Xi+Vi = e ±lC * r . ) are important (see Fig. 1), 
we represent the functions u(r) and v(r), similar to the 
one-dimensional case, as 



(P-Q)r- 



. (9) 



peFS, Px >o 



where Q = Q+ for wave vectors p v > and Q = Q_ for 
wave vectors p y < 0, respectively. 




FIG. 1: The Fermi surface 



where E g is the ground state energy and e n > is the 
energy of the n-th excitation. Following^ we obtain the 
eigenvalue equations 



*E U < 



For the doped case nesting vectors Q± are no longer 
equivalent. Therefore in the general case we consider 
vectors Q± as independent and make the substitution 



.(r + 5) + (^p(r)-fi)u a +m(r) au CT (r) "> 



Y A(r, r + 5; a)av a (r + S) 
s 



= e a u a (v), (7) 



- Y A*(r, r + 5; -a)au a (r + 5) + t ^ v a {r + 8) 

s 8 

-(yp(r) - tx)v a + m(v){-l) x '+y^v a (v) = e a v a (r), (8) 

where S = ±x, ±y. 

We suppose the d x i_ y i symmetry of the superconduct- 
ing order parameter A(r, r ± x; a) = erAd(r), A(r,r ± 
y; a) = — aA^r). The Fourier transform gives the usual 
dependence A p (r) = a J2$ A sc (r, r + 5; a) exp[— ip5] = 
2(cosp K — cospy)Ad(r). The system (jTj) - dHJ can be 
rewritten in the continuum approximation. Consider 
states near the Fermi surface (FS) (see Fig.l) and use lin- 
ear approximation for the quasiparticles spectrum. Since 
for SDW pairing components with wave vectors p and 
p — Q + (or p and p — Q_, where Q+ — Q_ = 27r(0, 1) 
is the lattice vector for the pure system without doping, 



m + (ri) exp(iQ + r + ) + m_(ri) exp(iQ_r_) + h.c. 

Eigenvalue equations ([7]), © take form similar to the 
ID case: = e*, with 



H 



-Ar, 



B r 



A_ r 





A-p + Q 



(10) 



774(1*) 



m±(r) 



(11) 

where 7?(r) = p - fp(r), <f T = (u p , u p _ Q , v p , v p -q) = 
(u + ,u-,v + ,v-), e p = -2t(cosp x + cosp y ) — fj,, V p = 
2t(smp x , sinpj,), and, as before, Q = Q + for wave vectors 
p y > and Q = Q for wave vectors p y < O.The sign in 
m± is taken by the same rule. 

For the case d x i_ y i symmetry we consider A_ p = 
A p = — A p _q = 2(cosp x — cosp y )A t ;(r), which corre- 
sponds to crA(r, r ± x; a) = — aA(r, r ± y; a) = A in the 
uniform ground state. We retained the main terms in the 
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expansion over A. In the higher order approximation, in- 
stead of the terms A_ p and A_ p+ q we would have to 
write A_p-i(VpA_ p )V r and A_ p+ Q-i(V p A_ p+ Q)V r . 
The continuum approximation is not valid for a band fill- 
ing very close to the half-filled case (the number of par- 
ticles per one site p — 1), where the Fermi velocity tends 
to zero at points p = (0, ±7r), (±7T, 0). 

In the homogeneous case p(x) — const, m, A^ = const 
for coexisting spin- and superconducting order parame- 
ters, so the eigenvalue spectrum has the form 



£ 2 = ( v /m2+ e 2(p)±^) 2 + A2, (12) 

with A p = Ad (cos — cos p y ). 

The self-consistent conditions are derived by substi- 
tution of functions u, v into (j4]), similar to the one- 
dimensional case. In the continuum approximation they 
read: 

p(r) = 2^[«u+ + u*_u.)f + (v* + v+ + v*_V-)(l - /)] 

(13) 

(_l)*«+» m ( r ) =4U[J2u*_u + f-J2vlv + (l-f)} (14) 

e e 

Aq(r) = 2g^2(v* + u + - v*_U-)[(l - f)(coa(p x - q x ) 

C 

+ cos(p y - q y )) - f{{cos(p x + q x ) + coa(p y + q y ))],(15) 

where / = l/(exp[e/T] + 1). We omitted spin indices 
since in our representation for wave functions all equa- 
tions are diagonal over spin. 



III. SPIN-CHARGE DENSITY WAVE 
STRUCTURES 

In the low doping limit the ground state of the model 
is the periodic charge-spin superstructure with the ab- 
sence of superconductivity: A = 0. Consider different 
structures, having close ground state energies. In real 
systems the exact ground state must be determined by 
taking into account real long-distance 3D interactions. 



A. Diagonal stripes 

For diagonal stripes we search the solution in the form 



where V p = 2t sinp x . The only difference from the consid- 
ered one-dimensional model^ is the dispersion of the ve- 
locity V p . This system is exactly solvable. In the ground 
state, at p — 1, we have m(r + ) = mo. Increased doping 
leads to the stripe structure. The one stripe solution has 
the form 



m(r_) 



mo tanh 



(17) 



where the width £ is defined from the minimum of the 
total energy. Solution (TIT)) corresponds to p = 1 in the 
thermodynamic limit for number of holes per lattice site. 
In our case (fl7j) is valid only in the vicinity of each single 
stripe that enters a periodic superstructure called stripe- 
phase (compare^ 9 -). Distinct from the Peierls model, 
where £ = Vp/mo, Vp — const, the present model has 
a more complicated spectrum. Besides continuum bands 
E 2 — V 2 k 2 + niQ we find some discrete levels (for a given 
p x ) inside the gap: 



E 2 n = ml\n(2 - An), 



(18) 



where n is integer number, < n < 1/A, and A = X(p x ) = 
V p /(£mo), or A = V p /V with £ = V/m . Each level in- 
side a gap forms a band due to dispersion of the coeffi- 
cient \{p x ). For A > 1 we obtain only one level, E = 0, 
with wave function ip± = (u + ± uJ)j\[2: 

1 

ip± oc 



(cosh r + /t;) l / x ' 

The wave functions of all states are described in terms 
of the hypergeometric function F(a,b\c\z), and for local 
levels they have polynomial form: tp±,n ~ 

(coshJo^-" ^ + 1 '"' ^ (1+tanh | )] 



For 1/2 < A < 1 two levels n = 0,1, can already ex- 
ist. Allowing for equality A = V p /V we conclude that 
both possibilities A < 1 and A > 1 take place, each one 
in the proper interval of p x . Similar to the ID case 14 
the appearance of the kink in the spin channel is accom- 
panied by the local charge distribution p(r+) — (p) ~ 
1/ cosh 2 (r + /£). An increase of the doping leads to the 
periodic spin-charge density superstructure. In the limit- 
ing case of "overdoping" (| p— 1 |>- l/£) the spin-charge 
structure becomes harmonic 

m(r+) oc sin(-7r | p— 1 | r + ), p(r + ) — (p) oc cos(27r | p—1 | r + ). 



up(r) = u p (r+), u p (r) = v p (r+), 

where r± — (x ± y)/ y2, p € FS. Substituting into Eqn. 
(|lip we obtain a one-dimensional eigenvalue equation 

d U 

- iVp- — u+ + —p(r + )u + + m{r+)u- = Eu+, 
or + 2 

d U 

m*(r+)u+ + iVp-^—u- + —p(r+)U- = Eu_, (16) 



B. Vertical stripes 

Using a similar ansatz for discrete equations 

u p (r) = u p (x), v p (r) = v p (x) 

we obtain a system of equations with V p — 2tsinp x , 
which is similar to the diagonal case. For the same values 
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mo, £, the parameter A in the considered case is less than 
for diagonal stripes. Therefore the condition n < 1/A 
can be valid for larger values of n, resulting in additional 
bands inside the gap, as it is seen from numerical results. 



C. Checkerboard structure 

As we have seen the spin-charge density structure may 
be arranged in vertical (horizontal) or diagonal direc- 
tions. Consider the solution with square symmetry. In 
the same approximation as before we find the solution of 
system (fTTj) in the form m±(r) = m(r±), u p (r) — u p (r±), 
v p (r) — v p (r±). Equations are decoupled and we obtain 



,„ du + . 

tv„— h m{r±)U- 

ar± 



m*(r±)u + + iV f 



du- 
dr± 



-- Eu^ 
Eu- 



in) 
(20) 



with V p — 2tsinp x , r± = (±x + y)/y2. The one "cross" 
solution has the form 

m + = mo tanh — , to_ = too tanh — , (21) 

The spectrum E and wave functions are found as above 
for the case of stripes. In the case of high doping the one 
kink solution is transformed to the periodic structure 

(S z (r)) cx (-ir+fm cos[7r(Vp- l)x] cos[tt(Vp- l)y], 

in which we considered the squared Fermi surface approx- 
imation with electron density p = |Q| 2 /27r 2 . 



IV. SUPERCONDUCTIVITY AND 
SPIN-CHARGE MODULATION 

A. Vortex solution 

Consider pure superconducting state (A(r) = 0). The 
BdG equations are decoupled. The first pair is 



iVpV r Up(r) + AptJp 



A*Mp + iVpV r Vp(r) = ev p . 



(22) 



(23) 



When the filling p is close to 1, the Fermi surface has 
nearly square form, therefore V p V r « V p d/dr±, depend- 
ing on signs p x , p y In this case the system of equations 
([22]) . (|23]) has the following vortex solution: 



A-p(r) 



sinh • 



i sinh t 



sinh 2 ^ + sinh 2 j- + 1 



(24) 



where A p = Ao(cos(p x ) — cos(p«)). For r_ = the 
order parameter has a kink form A p (r) cx tanhr + /£ s . 



For the case r + = the order parameter acquires 
the phase: A cx exp(j7r/2) tanhr_/£ s . In the di- 
agonal direction r+ = r_ the solution A p (r) cx 



tanh r+/£/y tanh 2 r + /£ s + 1 exp(i7r/4) has the phase 
7r/4. It is known that in one-dimensional case finite- 
band solutions of equations (|22|) - (|23|) are related to 
the soliton (kink) solutions of the nonlinear Schrodinger 
equation (NSE). Note, that along the curve sinhr_/£ s = 
a cosh r + /£ s the order parameter acquires the form of 
of a general kink solution of the NES: A p (r) ~ (ia + 
tanha;/£ s )/v'a 2 + 1 with the localized state in the gap 
with the energy E = A p a/y/a 2 + 1. 



B. Coexistence of spin-charge structure and 
superconductivity 

Consider solutions of equations (|11[) in the supercon- 
ducting region. 

By analogy with ID case^ we use ansatz: 

v± = j±u T 

which takes place in the uniform case. The term U p(r)/2 
in equations can be eliminated by the shift of wave func- 
tions u,v — > u, uexpi$, VpV^ = Up(r)/2. Considering 
e(p) — p, = on the Fermi surface we obtain for the case 

m (r) = \m(r)\e lv , A p (r) = \ A p (r)\e lVs , ip, ip s = const, 

the solution 7+ = ±ie l(v ~ V3) , 7_ = ±ie- l( - v+v °\ and the 



system (Ti"Tj) acquires the form 



A(r)w_ = Eu + 
A*(r)u + + iVpVu- = Eu^ 



(25) 
(26) 



with A(r) = (|m(r)|±i|A p |)e* ¥ ', and, as before, to = m±, 
depending on the sign of p y . Equations (|23|) , (|2"rJj) are ex- 
act provided that phases ip, p s are constant or slowly 
varying in space functions. We show that inhomogene- 
ity of the superconductor order parameter leads to the 
the origination of the antiferromagnetic order parame- 
ter. Consider a ID geometry case: u = u(r + ), where 
assumption of constant phases is valid. The solution of 
Eqs. (|2"5|) . ([21?]) describing the coexistence of supercon- 
ductivity and spin-charge density ordering, compatible 
with self-consistent equations, has the form of two bound 
solitons of the nonlinear Schrodinger equation, see, for 
example^ 

r _ cosh(2KX + a) + cosh(c 2 ± 2/3i)/\X\ 
1,2 ~ p cosh(2 ra + Cl ) + cosh(c 2 )/|A| ' [ ' 

where ±A(p) are positions of local levels inside the gap, 
k = y^A 2 — A 2 , and exp i(3 — A + in. For superconduct- 
ing and spin order parameters we obtain 

A sc = A p (l - rtanha(tanh(^ + |) - tanh(^ - |))), 

(28) 
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m = m rtanha(tanh(-i + -) - tanh(-± - -)) (29) 

where averaged over Fermi surface: T =< T p >=< 
A p /yjAj + ml > p , £ =< Vp/iApy/fptauha) > p , and 
we use the parametrization for the local level A. 

A 2 / A 2 \ 

A 2 = A2 p 2 ml + (30) 

A P + m o\ cosh 2 a J 

values of mo, A p = Ad(cos p x — cosp y ), £ and a dimen- 
sionless parameter a are defined by the self-consistent 
conditions (fl~3]) . ()14|) . The solution describes the spin- 
charge stripe in superconducting phase. 



charge-spin structures, superconductivity, and their co- 
existence. The origin of spin-charge periodic state (which 
is responsible for the pseudogap) is due to the existence 
of flat parallel segments of the Fermi surface (nesting) 
at low hole doping concentrations. Effects of commensu- 
rability lead to a pinning of stripe structure at rational 
filling points \p — 1| = m/n. As a result, there is an 
exponentially small (for large n) decrease in the total 
energy of the order dE exp(-constn) at any commen- 
surate point, stabilizing stripes, as in ID systems. For 
this reason, we think, stripes are mostly observable near 
n = 8 point (\p — 1| = 1/8). An increase of doping leads 
to the decrease of flat segments of the Fermi surface and 
attenuation of spin-charge structure. 



-Oh i 

4 Z g 4 N, 



FIG. 2: False colour plot of the coexisting superconducting 
(donward) and antiferromagnetic stripe-like (upward) orders. 
The envelope functions are plotted in real space, x and y 
coordinates are measured in units of correlation length £, a — 
5, A d = 1, m = 0.2. 

The spin inhomogeneity generates the charge distribu- 
tion Sp(r) ex. m 2 (r). Note, that two-soliton solution in the 
similar form was used for describing polaron-bipolaron 
states in the Peierls dielectrics^. 

The correlation length is increased in comparison to 
clean superconductor case as 

t = Uy/i + QJ&F (31) 

V. DISCUSSION 

We considered a simple self-consistent 2D model on 
a squared lattice to describe different states, including 



We found the solution describing the coexistence of 
superconductivity and stripes (28), (29). The decrease 
(or a deviation from the homogenous value ) of the su- 
perconducting order parameter generates the spin-charge 
periodic structure in this region. Note, that due to sym- 
metry of Eqs. (25)-(26) (duality A <-> im) we can write 
the same equation, describing the origin of superconduct- 
ing correlations in the region of a inhomogeneity of spin- 
charge density. The situation is qualitatively similar to 
the ID case^. Experimental data in underdoped high-T c 
cuprates LSCO 5 indicates that antiferromagnetic stripe- 
like spin-density order can be induced by magnetic field 
perpendicular to the CuO planes in the interval of fields 
much smaller than upper critical field H C 2 . The size of 
the magnetically ordered domains exceeds superconduct- 
ing vortex's core size £ s and the inter-vortex distance in 
the Abricosov's lattice. Our present theoretical results 
demonstrate that this is indeed possible in the simple 
Hubbard t-U-V model that we consider. In particular, 
the dimensionless parameter a in Eqs. (28)-(29) is an 
independent variational parameter and depends on the 
magnetic and superconducting coupling strengths^, as 
well as on the magnitude of the external magnetic field. 
Hence, the size ~ a x £ of the antiferromagnetic domain 
(see Fig. 2, upward red plane bump) can exceed the su- 
perconducting (and magnetic) Ginzburg-Landau correla- 
tion length £ when a(H) >> 1. Previously coexistence of 
superconducting order and slow antiferromagnetic fluctu- 
ations was studied merely on the basis of a phenomeno- 
logical Ginzburg-Landau free energy functional approach 
iniii. We note also, that equations Eqs. (25)-(26) can be 
simply extended to include d-density waves (DDW). 
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